clear all
close all

load resMat.mat


figure(1)

subplot(2,4,1), plot(tr,kr,'linewidth',2)
title('Capital k')
subplot(2,4,2), plot(tr,Hr)
title(texlabel('Cumulative knowledge H'))
subplot(2,4,3), plot(tr,etar)
title('Shadow price of H eta')
subplot(2,4,4), plot(tr,lambdar)
title('Goods shadow price lambda')
subplot(2,4,5), plot(tr,invr)
title('Capital investment i')
subplot(2,4,6), plot(tr,jr)
title('Technology investment j')
subplot(2,4,7), plot(tr,pr)
title('Renewable energy cost p')
subplot(2,4,8), plot(tr,epsrelr)
title('Energy relative price epsilon/lambda')

figure(2)


subplot(3,3,1), plot(tf,kf)
title('Capital k')
subplot(3,3,2), plot(tf,Sf)
title('Fuel use S')
subplot(3,3,3), plot(tf,Nf)
title('Fossil technology capital N')
subplot(3,3,4), plot(tf,lambdaf)
title('Goods shadow price \lambda')
subplot(3,3,5), plot(tf,sigmaf)
title('Fossil fuel user cost \sigma')
subplot(3,3,6), plot(tf,nf)
title('Fossil technology investment n')
subplot(3,3,7), plot(tf,invf)
title('Capital investment i')
subplot(3,3,8), plot(tf,g)
title('fossil energy cost g')
subplot(3,3,9), plot(tf,epsrelf)
title('Energy relative price \epsilon/\lambda')
% 
% yr = A*kr;
% yf = A*kf;
% y = [yr;yf];
% 
% costf = g.* yf;
% costr = pr.*yr;
year = [tr;tf];
conShare = [cr;cf]./y;
epsrel = [epsrelr;epsrelf];
invShare = [(invr+jr);(invf+nf)]./y;
figure(3)

subplot(1,3,1), plot(year,conShare)
title('Consumption Share of Output c/Ak')
subplot(1,3,2), plot(year,epsrel)
title('Energy Relative Price eps/lambda')
subplot(1,3,3), 
plot(year,invShare)
title('Total Investment Share of Output (i+n&j)/Ak')
